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Recent experiments on the quasi 1-D antiferromagnet Cu 
Benzoate revealed a magentic field induced gap coexisting 
with (ferro) magnetic order. A theory explaining these find- 
ings has been proposed by Oshikawa and Affleck. In the 
present work we discuss consequences of this theory for inelas- 
tic neutron scattering experiments by calculating the dynam- 
ical magnetic susceptibilities close to the antiferromagnetic 
wave vector by the formfactor method. 
PACS: 74.65. +n, 75.10. Jm, 75.25,+z 



I. INTRODUCTION 

It has been known for some time that Cu Benzoate 
is a quasi-lD S = 1/2 (Heisenberg) antiferromagnet [Q. 
In a recent neutron scattering experiment Q in a mag- 
netic field, the existence of field-dependent incommen- 
surate low energy modes was established. The incom- 
mensurability was found to be consistent with the one 
predicted by the exact solution of the Heisenberg model 
in a magnetic field. However, the system exhibited an 
unexpected excitation gap induced by the applied field. 
These findings were quite surprising as the coexistence of 
a gap and a magnetization is inconsistent with rotational 
symmetry around the direction of the magnetization || . 

The resolution of this puzzle goes as follows: Dender et 
al [|| suggested that the gap could be due to the staggered 
magnetic field generated by the alternating g-tensor in 
Cu Benzoate; Oshikawa and Affleck then showed that 
this mechanism (in particular the generation of a stag- 
gerd field perpendicular to the direction of the applied 
magnetic field), together with the Dzyaloshinskii-Moriya 
interaction present in Cu Benzoate, indeed can account 
for the experimental findings. They proposed that Cu 
Benzoate is described by the Hamiltonian 



(1) 



where H 3> h. The low energy effective theory of ([!]) is 
obtained by abelian bosonization and is given by a Sine- 
Gordon model with Lagrangian density pi 



C = -{d^f + \cos{(i 6) . 



(2) 



Here 9 is the dual field and fulfils d x § = —d t @, <9 f $ = 
-d x S. This implies that (d^) 2 = (<9 M 9) 2 and (§) thus 
is a Sine-Gordon model for the dual field. We note that 
all of our analysis is based on the model (|l|) and (||) so 
that the results carry over to any material described by 
this theory (with appropriate fit of -j to experiment). 



The coupling constant (3 is a function of the applied 
field H and is related to the "dressed charge" (see || 
and references therein) of the isotropic Heisenberg chain 
in a magnetic field by 



For further convenience we define the quantity 

5 8tt - /3 2 



(3) 



(4) 



The dressed charge can be determined exactly from the 
Bethe Ansatz solution. It is given by Z = Z(A), where 
Z(X) is a solution of the integral equation [H 



+ (a-m) 5 



Z{p) = 1 . (5) 



The integration boundary A is a function of the magnetic 
field H and is determined by the condition e(A) = 0, 
where 



£ (A) 



1 



df i 



2tt J_ a 1 + (A — (j,)- 



2H 



J A2 + i 



(6) 



These integral equations can easily be solved numerically. 
In Fig. |l| we plot j3 2 /2ir of (where/3 is the coupling con- 
stant of the Sine-Gordon model (||)) as a function of the 
applied magnetic field H for Cu Benzoate. We note that 
at zero field /3 2 /2n = 1. 
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FIG. 1. Coupling constant !3 2 /2tt as a function of the ap- 
plied magnetic field. 
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The curve in Fig. [I] terminates at a value of H sa 27 
Tesla (H = 2J) at which the underlying Heisenberg 
model experiences a transition to the saturated ferromag- 
netic state. 

The spectrum of the Sine-Gordon theory (j^) is well 
known || (see also [@||): it consists of a soliton- 
antisoliton doublet of mass M and their bound states 
which are called "breathers" . For Cu Benzoate the soli- 
ton mass as a function of the applied magnetic field was 
determined by Oshikawa and Affleck 



(7) 



where the induced staggered field h is determined as a 
function of the applied field H (and the Dzyaloshinskii- 
Moriya interaction in Cu Benzoate) in Q . 

In addition there are [|] (here [x] denotes the largest 
integer smaller than x) breathers, denoted by B%, B2, ■ ■ ., 
with masses 



M n = 2M sinmr£/2 ; n = 1, . . . , [l/£] 



(8) 



The mass spectrum as a function of magnetic field for 
CuBenzoate is plotted in Fig. || 
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FIG. 2. Breather masses in units of the soliton mass as 
functions of the applied magnetic field. 

We see that for small fields there are three breathers, 
one of which is just below the soliton-antisoliton contin- 
uum starting at 2M. For higher fields further breathers 
split off this continuum. However they appear at ener- 
gies higher than the threshold of the B\B\ two breather 
continuum. For small fields the dependence of the mass 



ratios 



M„ 
M 



on H is 



M„ ~ 2M 



sin nir/6 + 



9 In h/ Jh 



cos nir/6 



(9) 



The Lagrangian (g) can be brought to standard Sine- 
Gordon form by performing the duality transformation 
— > $. The bosonized expressions of the spin operators 
change under this duality transformation as well. We 
find [pFj|-p^[ S n — > aoS(x), x = nao, 

S(x) = J(x) + (-l) n n{x) , 

J z = \d x @ , 
p 

J+ = — !— exp cos f ^9 - 2Sx] , 

2na V P / 

n x (x) = cos(/3$(x)), 
na 

n y (x) = -^-sin(/9$(x)), 
na 

n*( x ) = — cos(^e(a;) - 2Sx). (10) 
nao p 

Here n a {x) are the components of the staggered magne- 
tization, J ±:Z are the current operators, nao = x and A 
is a nonuniversal coefficient. The quantity 5 is the differ- 
ence between the Fermi momenta for up spins and down 
spins. It is easily determined from the exact solution of 
the XXX Heisenberg chain in a magnetic field (see || and 
references therein) but is not important for the present 
work. 

The dynamical magnetic susceptibilities for wavevec- 
tors close to the antiferromagnetic wave vector n (it ± 2<5 
for the longitudinal part \ zz ) are then given by 

/DO />OC 
dx / dt e *(«+fc)*-i(«-T)* x 
-00 Jo 
x ([cos /3$(t, x), cos /?$(0, 0)]) 

/oc /*oo 
dx / dt e i(<-+fc)*-*(9-T)* x 
-oo Jo 
x([sin/3$(t,ir),sm/3$(0,0)]) 

/•GO /*00 

X zz (w,£7)cx / dx dt e <(«-He)t-ifo-ir)x x 



./-oo ./0 

x ([cos(^e(t, x) - 2(5x), cos ^6(0, 0)]). (11) 

The transverse susceptibilities \ xx an d X yy can be 
determined exactly through the formfactor approach to 
quantum correlation functions Jl3|-p^[ (our discussion fol- 
lows ]l6|). The longitudinal susceptibility involves the 
dual field which is a nonlocal operator and cannot be 
treated in the same way. We therefore restrict our analy- 
sis to the transverse susceptibilities. As the longitudinal 
modes become soft at the incommensurate wavevectors 
it ±25 rather than at the Neel wave vector it (where the 
transverse modes are soft) it is possible to separate \ zz 
from x aa j a — x,y experimentally by performing con- 
stant wavevector scans around q — tt ± 2<5 and q = n 
respectively. 



where ho = -\/87r 3 /e. 
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Insertin g a resolution of the identity in the correlation 
functions (|llj) we obtain formfactor expansions of the 
form 



oo „ 



dgi . . . d6>.„ 
(2n) n n\ 



(6i...9 n )f 



n— 1 

5[(7r-g)-X;^i8inh( 



cj — £^ • Mj cosh 9j + ie 
6[(n-g)+E j M j smh9 j ] 

U + Si -^3 C0Sn + * e 



(12) 



and an analogous expansion for x vy '■ The quantities F cos 
in ([l2|) are the formfactors of the operator cos/3<I> and 
have been calculated exactly [^3| - [l5| . The expansion jl^ ) 
is in terms of multiparticle states. In order to make con- 
tact with the neutron scattering experiments we are in- 
terested in the imaginary parts of the dynamical suscepti- 
bilities, which are proportional to the neutron scattering 
cross sections. Because all elementary excitations in the 
Sine Gordon theory have a gap, n-particle states will con- 
tribute to Smx""^! q) i a — y) only at energies larger 
than the n-particle gap. In other words, by taking into 
account only the first few terms in the expansion ( |l2| ) we 
obtain an exact expression for 3mx aa (u,g) for energies 
below a threshold that depends on the omitted terms. 

As the Sine-Gordon theory is Lorentz invariant the fi- 
nal answer for the susceptibilities cannot depend on lu 
and q independently but only on the combination 



s = y/uj 2 - (tt - q) 2 



(13) 



The formfactor expansion is now readily performed us- 
ing the results of |l3H15|l and the bootstrap procedure 
to determine breather formfactors. Using that cos (3$ is 
even under charge conjugation Jl4],|l6| we find 



[i/€] 



<3m X xx {uj, q) cx 2vr ^ Z 2„ S(s 2 - M. 



2 ) 

In) 



-Re 



\F cos [9(Mx, Mi, s^up 



+25ie 



s^/s 2 - AM 2 
\F cos [9{M,M,s)} + -\ 2 



where 



9(mi, mi-, s) — arccosh 



,Vs 2 - 4M 2 



s — mr — mi 



(14) 



(15) 



Here the dependence of £ on the applied field is given 
by (||) and Fig. |^. The first terms in (jl4|) correspond 
to single-particle breather states. The (squares of the) 
breather formfactors are given by 



Z-> 



2(sin2<) 2 
cot 7r£ 

exp 



) dx (sinh 2£a;) 2 sinh a; (1 — £) 



Z 4 = 



o x cosh x sinh 2x sinh £x 
2(sin4<) 2 



x 



(cot7r?) 2 (cot3</2) 2 

dx (sinh 4£ir) 2 sinh x(l — £) 
x 



exp 



-2 



j dj cosh a; sinh2x sinh £2; 
The soliton-antisoliton formfactor is 

) 2 cosh#/£ + cos7r/£ 



(16) 



\F+^\P)\ = 



x exp 



(2cot7r£/2 sinh( 



£ 2 cosh2(9/£- cos2tt/£ 

dx [cosh2a;cos2x0/7r — 1] sinhx(l — £) 



coshx sinh 2x sinh^a; 



(17) 



and the -Bi-Bi breather-breather formfactor is 
(sinh (9) 2 



X 



x exp 



where 



(sinhfl) 2 + (sinO 2 
dx cos 2x9 /tt sinh £x sinh x sinh x(l + £) 



(sinh 2x) 2 



A = 2 cos tt C/2V2sinX/2 exp 



^ dx x 



2tt sin x 



(18) 



(19) 



The next most important contribution (i.e. the term with 
the lowest threshold) in the expansion (|l4] ) comes from 
B1B3 breather-breather states. It will contribute at en- 
ergies larger than Mi + M3, where Mi 3 are given by 

©• a 

In Fig. [5| we plot our result for x^i 1 ^, q) as a function 
of s/M = yju 2 - (tt - q) 2 /M (where M is given by (0)) 
for an applied field of H = 3.5 Tesla, which is one of the 
values studied experimentally in [0. 
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FIG. 3. Transverse susceptibility x* x ( s ) f° r H = 3.5Tesla. 
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We see that there is one coherent mode at s w 1.5M 
with weight Z 2 ~ 1.02. At s w 1.68M the B1.B1 con- 
tinuum appears and at s = 2M the soliton-antisoliton 
continuum. The susceptibility is regular at both thresh- 
olds. 

Fig. |] shows x xx {u, q) as for H = 7 Tesla. The picture 
is qualitatively similar to the one at 3.5 Tesla although 
the soliton-antisoliton threshold has become more pro- 
nounced. The weight of the delta-function correspond- 
ing to the breather B 2 has diminished to Z' 2 ~ 0.97. As 
the field increases a singularity develops at the s = 2M 
threshold until at H ss 9.6 Tesla the breather B4 splits off 
the continuum. The situation at H = 14 Tesla is shown 
in Fig. | The weight of B 2 is Z'J, a 0.85 whereas the 
fourth breather has only a very small spectral weight of 
Z 4 « 0.01. 
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FIG. 4. Transverse susceptibility % XX { S ) f° r H = 7Tesla. 
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FIG. 5. Transverse susceptibility x**^) f° r H = 14Tesla. 
Let us now turn to Sm^'fw,?). It is clear from the 



form of the Hamiltonian ([!]) that this will be different 
from Qmx xx (w, q). The formfactor expansion yields 



[i/«] 



%m X vv (u>, (?) oc 2tt £ Z 2n-x S(s 2 - M| n _ 1 ) 



n=l 



+25Re 
+23f?c 



|F sm [6»(M, M,s)]+_| 



4M 2 



|F sm [#(Mi, M 2 , s)] 



12 



yf{s 2 - M'{ - Af|) 2 - iMfAq 



(20) 



Here the breather formfactors are given by 
8(costt£/2) 4 



Zi 



cot tt£/2 
exp 



^ , dx sinh£a; sinhx(l — £) 



/ a; cosh a; sinh 2a; 
4sin3< (sin37r^/2) 2 



exp 



(cot<) 2 

t , dx (sinh3£a;) 2 sinha;(l — £) 



/q x cosh a; sinh 2a; sinh£a; 
The soliton-antisoliton formfactor is found to be 

(2cot?r£/2 sinh6») 2 cosh 0/£ - cos ?r/£ 



(21) 



x exp 



£ 2 cosh 26>/£- cos 2tt/£ 

da; [cosh2a;cos2a;6'/7r — 1] sinh a; (1 — £) 



cosh a; sinh 2a; sinh £3; 



(22) 



Finally we take into account the B\B 2 breather-breather 
state which has a formfactor of 



X ( 1 + 



tan 7r£ 



.g(0-wr£/2).g(0 + M r£/2)| 2 A 6 



x exp 
x exp 

where 



4 cos 7r£/2(cosh6> + cos 7r£/2) 

dx cos 2x8 /it sinh 2£a; sinh a; sinh a; ( 1 + £) 



g a; (sinh 2a;) 2 

cfa; cosh 2a;£ sinh £3; sinh x sinh x(l + £) 
a; (sinh 2a;) 2 



9( x ) = 



sinh : 



sinh x — i sin 7r£ 



(23) 



(24) 



The next most important contribution to (|20j ) is due to 
B1B1B1 three breather states with a threshold at 3Mx. 

In Figs. ^[]| we plot x yy ( s ) f° r three different values 
of the applied field, namely H = 3.5, 7, 14 Tesla. The 
emerging picture remains qualitatively unchanged for all 
three values of H. There are two coherent modes corre- 
sponding to the breathers B\ and B3. At s = 2M the 
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soliton-antisoliton continuum starts and at s = Mi + M2 
we observe the onset of the B\Bi breather- breather con- 
tinuum. The latter one exhibits a singularity as a func- 
tion of s. The weight of the delta functions correspond- 
ing to the breather B\ decrease with increasing field: 
Zi ps 2.07, Z[ ps 2.04 and Z'{ = 1.95. For B 3 we find 
a decrease with increasing field after an initial increase: 
Z 3 PS 0.34, Z' 3 ps 0.37, Z% = 0.36. Note that the spectral 
weight of the heavier B3 particle is always significantly 
smaller. 



1.0 



3 



0.0 



0.0 0.5 1.0 1.5 2.0 2.5 

s/M 



FIG. 6. Transverse susceptibility x vv ( s ) f° r H = 3.5Tesla. 
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FIG. 7. Transverse susceptibility x vv ( s ) f° r H = 7Tesla. 
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FIG. 8. Transverse susceptibility x vv (s) for H = 14Tesla. 

Let us now compare our results with the experimen- 
tal findings. Our results for the imaginary part of the 
dynamical susceptibilities are consistent with the exper- 
iment at 7 Tesla (Fig. 3 (c) of ||): as the beam was 
unpolarized the experiment observed \{x xx + X vv ) (the 
contribution from \ zz emerges only at energies greater 
than gfiBH). The observed peaks at 0.17meV, 0.34meV 
and 0.44meV correspond to the three breathers B\, B2 
and £?3 and agree very well with our prediction for the 
mass spectrum (||) and (^|) (the first two peaks were al- 
ready discussed in B). In order to compare our calcu- 
lated weights of the delta-functions in the susceptibilities 
to the experimentally observed intensities one needs to 
convolve with the instrumental resolution. Convolving 
our exact results with a Gaussian we obtain a fit to the 
experimental data shown in Fig. |^. 
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FIG. 9. Calculated (solid line) vs measured (data points 
taken from Ref. [2]) intensity at H = 7Tesla. 

Wc sec that the theory is in good agreement with ex- 
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periment. 
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